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Abstract 

Atmospheric wind speeds and their fluctuations at different locations (on- 
shore and offshore) are examined. One of the most striking features is the 
marked intermittency of probability density functions (PDF) of velocity 
differences - no matter what location is considered. The shape of these 
PDFs is found to be robust over a wide range of scales which seems to 
contradict the mathematical concept of stability where a Gaussian distri- 
bution should be the limiting one. 

Motivated by the instationarity of atmospheric winds it is shown that the 
intermittent distributions can be understood as a superposition of differ- 
ent subsets of isotropic turbulence. Thus we suggest a simple stochastic 
model to reproduce the measured statistics of wind speed fluctuations. 

1 Introduction 

Atmospheric wind may be seen as a prime example of a turbulent velocity field 
with very high Reynolds numbers of about Re ~ 10^ Reynolds numbers as 
large as this prevent analytical calculations and direct numerical simulations. 
Therefore the flow has to be described in a statistical way. For the estimation 
of extreme loads as well as for risk estimations the statistics of velocity fluctu- 
ations u{t) and velocity differences should be known. It has been shown that 
these statistics obey non-Gaussian, intermittent distributions (e.g. [S]) that 
directly correspond to an increased number of wind gusts ^ . 
Nevertheless, for most technical and meteorological problems fluctuations as 
well as fluctuation differences are assumed to obey Gaussian statistics. There- 
fore simulations of atmospheric velocities are often based on Gaussian processes 
0. 

Fluctuation differences are commonly measured using velocity increments: 

Urit) := u{t + t) - u{t) . (1) 

Large increment values can be identified as wind gusts as long as the time step 
r is rather small. The demand for a small r-value (typically less than a minute) 
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is due to the fact that gusts are related to large velocity rises during short times. 
Rises that occur over time steps of several hours - for instance - are not called 
wind gusts but large scale variations. 

The difficulty to fix a suitable time scale mirrors the fact that atmospheric winds 
exhibit variations on any time scale - in principle ranging from seconds (and 
less) up to centuries. For most practical applications, such as engineering and 
meteorology one mainly distinguishes between large scale variations such as di- 
urnal, weekly and seasonal changes and variations on small scales often referred 
to as atmospheric turbulence or gustiness [J. The existence of a mesoscale gap 
as proposed by which divides small (micro) and large (macro) scales in a 
more rigorous way has strongly been debated in recent years (e.g. 0, [HJ). 

In this paper we focus on the scale dependent statistics of atmospheric in- 
crements and compare them to that of homogenous, isotropic and stationary 
turbulence^ as realized in laboratory experiments. For isotropic turbulence the 
statistical moments of increments, the so-called structure functions have been 
intensively studied ^U]. Their functional dependence on the scale r is described 
by a variety of multifractal models. Besides the analysis of moments, probability 
density functions (PDFs) are often considered. These show a transition from 
Gaussian distributions to intermittent (heavy-tailed) ones as scale decreases. 
Unfortunately the analysis of moments as well as that of probability density 
functions p{ur) is more or less restricted to the inertial range - the range of 
scales larger than the Taylor scale Q (where dissipation effects become signifi- 
cant) and smaller than the integral scale^ T. 

The challenge is to describe and to explain the measured fat-tailed distributions 
and the corresponding non-convergence to Gaussian statistics. Large increment 
values in the tails directly correspond to an increased probability^ (risk) to ob- 
serve large and very large events (gusts). As pointed out in [T^ the probability 
to observe large events - e.g. events twice as large as a reference - can become 
negligible for a Gaussian while for heavy-tailed distributions there is still a sig- 
nificant probability to observe it. 

The atmospheric PDFs - we examine here - differ from those of turbulent lab- 
oratory flows where - with decreasing scale - a change of shape of the PDFs 
is observed (e.g. 0). For large scales the distributions are Gaussian while for 
small scales they are found to be intermittent. The atmospheric PDFs how- 
ever change their shape only for the smallest and then stay intermittent for 
a broad range of scales. Such a constant shape for larger and larger scales is 
expected only for stable distributions such as Gaussian ones or the Levy stable 
laws Although the decay of the tails indicates that distributions should 

approach Gaussian ones (as for isotropic turbulence) they show a rather robust 
exponential-like decay. This point will be clarifled in chapter 13.21 

^For simplicity we will use the term 'isotropic turbulence' instead of homogenous, isotropic 
and stationary turbulence. 

^Normally Q and T denote length scales. For constant mean velocities and applying Tay- 
lor's hypothesis of frozen turbulence length- can be defined as time-scales as well. Here we 
will proceed with corresponding time scales. 

^The probability to observe an increment Ut e [ut,Ut + dur] is just given by p(ur)duT- 
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In chapter b.ll and lT^ it will be shown that atmospheric increments behave quite 
similar to those of isotropic turbulence for small scales but differ significantly 
for large ones. We therefore introduce a model - chapter 1^.31 - that interprets 
atmospheric increment statistics as a large scale mixture of subsets of isotropic 
statistics. When mixing is weak the same statistics as for isotropic turbulence is 
recovered while for strong mixing robust intermittency is obtained. In chapter 
13] the results are briefly discussed. 



The analysis presented in the following is based on one laboratory and four 
different atmospheric data sets. In addition to accessibility reasons the latter 
were chosen in such a way that their environmental and meteorological charac- 
teristics differ significantly. Therefore one offshore and three onshore data sets 
are examined. Additionally a laboratory data set was chosen as an example of 
an approximately isotropic turbulent wake-flow. 

The first data set was recorded in October 1997 near the German coastline 
of the North Sea in Emden at a height of 20 m by means of an ultrasonic 
anemometer ^Sl- The wind speed was measured continuously over a period of 
275 hours. In the following this data set will be referred to as Onl. 
The second data set - denoted as On2 - was obtained from a hot-wire mea- 
surement 6 m above the ground (in flat terrain) Here the wind speeds of 
approximately one hour were considered. 

OnS - the third data set was recorded in a very complex terrain near Oberzeir- 
ing (1900 m above the sea level) in Austria in 2001 J^l- The velocity was 
measured by means of an unltrasonic anemometer. The data consists of 255 
non-successive blocks of 4 hours length. The choice was made in order to obtain 
complete and continuous data within each block. 

The fourth data set - referred to as Off - was recorded during an offshore mea- 
suring campaign at Roedsand in the Danish Baltic Sea at 30 m height between 
1998 and 1999 ^16 . From this period 58 non-successive days were chosen. Again 
the choice was made in order to obtain complete and continuous data for each 
day. 

The laboratory data - denoted as Lab - were obtained from a wake flow mea- 
surement in the wind tunnel of Erlangen in 1998 |17j . A hot-wire was located 
2 m behind a cylinder of diameter D = 0.02 m in the plane of the cylin- 
der. Here the turbulent flow can be considered to be locally isotropic and fully 
developed. With a mean velocity of n = 20.9 ms~^ a Reynolds number of 
Re = uDu^^ ~ 30, 000 is obtained. Taylor and integral scale are found to be 
~ 2 • 10~^ s and T ~ 6 • 10~^ s respectively. 

For atmospheric data sets it is difficult to define an integral scale T because of 
the instationarity of atmospheric velocities that causes very long-range correla- 
tions R{t) oc {u{t + r)n(t)) so that the integral time 
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cannot be estimated properly. Only for On2 an estimate of the integral time 
can be given because of the rather constant flow conditions during the short 
measuring period of about 1 hour. 

For On2 the Taylor scale Q can be calculated while for all other atmospheric 
data sets the sample frequency is too low for 6 to be resolved. 
In Tabled a short overview of the most important specifications of all five data 
sets is given. 





Onl 


On2 


On3 


Off 


Lab 


N[l] 


3,958,874 


20,480,000 


14,688,000 


25,056,000 


12,500,000 


blocks [1] 


1 


1 


255 


58 


100 


/ [Hz] 


4 


5,000 


4 


5 


100,000 


u [ms ^] 


3.4 


8.3 


6.6 


9.6 


20.9 


min [ms^^] 


0.0 


0.8 


0.02 


0.6 


15.9 


max [ms^^] 


18.1 


18.1 


39.0 


36.1 


25.5 


W [ms^^] 


1.7 


2.3 


4.2 


3.2 


1.1 


e[s] 


< 0.25 


w 0.01 


< 0.25 


< 0.2 


2 • 10""^ 


T[s] 




« 14 






« 6 • 10"^ 



Table 1: The table summarizes some characteristic values of the different data sets, from top 
to bottom these are number of data points N , number of blocks, sample frequency f , mean 
velocity u, minimum, maximum, average standard deviation ct, Taylor scale Q and integral 
scale T. 

3 Analysis 

3.1 Scaling in Isotropic and Atmospheric Turbulence 

The central assumption for turbulent velocity time series is that they have a self- 
similar structure (a direct consequence of scale invariance of the Navier Stokes 
Equation) in the inertial range. This means that within this range the disorder 
of velocity fluctuations has a similar structure on every scale but with a scale- 
dependent magnitude. To quantify this one usually calculates the (absolute) 
moments of increments 

^? = (Kn , (3) 

which are also called structure functions of order n. Instead of calculating 
the absolute moments one can also consider (n") (see jlB, for a more detailed 
discussion). 

In isotropic turbulence structure functions are assumed to scale as: 

5;'ocr^" . (4) 

A linear (monofractal) scaling exponent corresponds to a self-similar struc- 
ture as proposed by Kolmogorov in 1941 jl9j who found that the scaling expo- 
nent should be 
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due to dimensional reasons. Instead of this linear behaviour various experi- 
ments suggest that the scaling exponent is a non-linear function of n. In 1962 
Kolmogorov j2()j introduced the following non-linear (multifractal) exponent 

Cn = ^ + ^(3n-n2) (6) 

motivated by the model of a turbulent cascade with a log-normally distributed 
energy transfer rate. The parameter // is called intermittency correction and is 
found to be close to 0.25 [2J which corresponds to = f — M = 1-75. This is 
in agreement with the examined data sets where ^6 takes values between 1.67 
and 1.80 as shown in Fig. ^). 



Figure 1: Scaling Exponents and Structure Functions 




a) This plot shows the scaling exponents 
(n as a function of order n for all five 
data sets. Additionally the linear scaling 
law given by Eq. (straight line) and 
the non-linear laws according to Eq. Q 
(curved line) and Eq. Q (curved dashed 
line) are shown. 



b) The structure functions of order 2, 4 
and 6 are plotted against that of order 
3 in a double- logarithmic presentation 
(vertically shifted for clarity of presenta- 
tion). The open symbols belong to Onl, 
the filled ones to Lab. 



In [22j another formula was proposed which seems to fit experimental data 
slightly more accurately than Eq. ©: 



There are other models besides these two multifractal ones, e.g. \2'A\ . |24j . The 
differences in ^„ in all these models are rather small (at least for small orders) 
so for simplicity we will restrict following discussions to the models given by 
Eq. © and Eq. ©. 

To estimate the dependence of Cn on n one has to calculate C„ first. The most 
common way to do this is to plot log{SJ^) against log{S^^) - a method referred 
to as Extended Self Similarity (ESS) 25 . The slope of the resulting line is 
equal to Cn- This is shown in Fig. ^) exemplary for data sets Onl and Lab. 
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Figure 2: The open symbols show the flatness F of Onl as a function of scale with r = 0.25 s. 
The filled symbols represent the flatness of Lab with r = 0.002 T ~ 10~^ s. Additionally fits 
according to Eq. ^ (straight lines) are drawn in with fi = 0.29 (Onl) and /i — 0.26 (Lab). 
The horizontal dashed line marks F = 3 (as for a Gaussian distribution). 



In both cases the slopes are in quite good agreement with Eq. ^ and Eq. (O. 
The slopes of Onl show small deviations from linear behaviour for large values. 
These correspond to large scales r that might not belong to the inertial range 
anymore. This already indicates that care should be taken when transferring 
standard analysis of isotropic to atmospheric turbulence. 

The difference between isotropic and atmospheric turbulence statistics becomes 
more obvious when calculating the flatness F. Assuming inertial range scaling 
according to Eq. Q and the flatness should scale as well and is given by: 

F := oc rC-^C^ . (8) 

If Eq. ((HI) is a suitable description^ the flatness scales according to: 

F oc T-^'^/^ , (9) 

as is easily shown inserting Eq. ^ into Eq. (jH)). As shown in Fig. |21 the 

measured flatness of the data shows this scaling behaviour but the absolute 
values of F are very different for different data sets. While for the Lab data set 
flatness approaches F « 3 for large r it saturates at F ~ 6 for the Onl data set. 
The flatness of OnS and Ojff saturates at ~ 5 while for On2 it goes down to 3.5. 



Calculating the flatness of a variable x is often done to estimate the shape of 
the PDF p{x). A Gaussian distribution has flatness 3. Deviations from this 
value can be taken as a hint for a non-Gaussian shape of PDFs. In the following 
the increment PDFs of the given data sets will be examined. 

*Eq. Q and other multifractal models yield very similar results because such low-order 
exponents as (4, and ("2 are quite indistinguishable from each other. 
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3.2 Probability Density Functions in Turbulence 



Alternatively to the analysis of scaling exponents one can directly investigate 
the PDFs of velocity increments p{ur)- The scaling behaviour as well as all 
moments (including derived quantities such as flatness or skewness) are imme- 
diately given if the distributions on every scale are known. 
In principle the knowledge of all moments and the knowledge of the PDFs 
should contain the same information as can be seen by means of the charac- 
teristic function ip{k) - which is the Fourier transform of the PDF - defined 



as: 



n=l 



(10) 



Nevertheless the relation between moments and PDFs is not unique. In it 
is pointed out that two different PDFs can have exactly the same moments. 
From many experiments of isotropic turbulence it is well known that the shape 
of a PDF changes with scale. Going from larger to smaller scales the dis- 
tributions become more and more heavy-tailed while for r > T a Gaussian 
distribution is obtained. This scale-dependent shape corresponds to non-linear 
scaling exponents - as introduced in Eq. (jB)) and Eq. ((T)) - while a linear be- 
haviour 5" = ariT"" = a^/?'^ according to Eq. Q leads to a constant shape. 
This can be seen by means of the characteristic function in Eq. that stays 
the same for a linear exponent only k is rescaled according to k = f3k. 
The analysis of scaling exponents thus focusses on the change in shape of distri- 
butions while the shape itself is of minor interest and could even be determined 
wrongly as shown at the end of the last chapter. Therefore we will henceforth 
focus on the analysis of PDFs. 

In accordance with Eq. ® (that was derived from the assumption of a log- 
normally distributed energy transfer rate) B. Castaing et al. |9j introduced a 
model in which the increment distribution p{ur) is interpreted as a superposition 
of Gaussian ones p{u-r\(j) with standard deviation a. The standard deviation 
itself is distributed according to a log-normal distribution /(cr). The increment 
distribution thus reads 



P{Ut 



j da p{ur\a) ■ f{a) 



da 



1 



a^J2^T 



exp 



K 

2(j2 



1 



aXV^ir 



exp 



ln^{a/ao) 
2A2 



(11) 



and will henceforth be referred to as the Castaing distribution. 
Two parameters enter this formula, namely ctq and A^. The first is the median 
of the log-normal distribution, the second its variance. The latter determines 
the form (shape) of the resulting distribution p{ur) and is therefore called form 
parameter. On one side the larger becomes the broader the log-normal 
distribution and the broader the range of a that contributes to the integral in 
Eq. 1)11(1 . On the other side the range of a becomes smaller and smaller with 
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Figure 3: 

P(ut) [a.u.] 



10' 




-5 5 



a) Symbols represent normalized PDFs 
(with their scale-dependent standard 
deviation = of the Lab data 

set. From top to bottom t takes the val- 
ues: 0.01 T, 0.05 T, 0.2 T and 1.0 T 
where T ~ 0.005 s denotes the inte- 
gral time. Straight lines correspond to a 
fit of the distributions according to Eq. 
GUI- 



Density Distributions 



P(Ut) [a.u.] 




b) Symbols represent normalized PDFs 
of the atmospheric data set Onl with 
r = 0.5 s, 2.5 s, 25 s, 250 s and 4000 s. 
Straight lines correspond to a fit of the 
distributions according to Eq. I|ll|l . All 
graphs - in (a) and (b) - are plotted 
in a semi-logarithmic presentation and 
are shifted against each other in vertical 
direction for clarity of presentation. 



decreasing form parameter. In the limit of vanishing A the log-normal becomes 
a delta distribution 

= S{a-ao) , (12) 

so that p(ur) is reduced to a Gaussian distribution with variance ctq. 
With a proper choice of the form parameter the PDFs p{ur) can well be fitted 
as it is shown in Fig. OK) and Fig. Wp). In Fig. OK) - where the Lah data 
set is presented - the expected change of shape from intermittent to Gaussian 
distributions with increasing scale is clearly seen. 

In contrast to this behaviour the PDFs of Onl in Fig. ISb) look totally different. 
They are much more intermittent and do not approach a Gaussian distribution 
even for very large scales. For scales larger than about 25 s the shape remains 
rather constant. This is in accordance with the finding of the slow decrease of 
flatness shown in Fig. |21 because in the log-normal model flatness is linked to 
the form parameter j27j according to 

A^ oc In 

In this sense constant flatness larger than 3 corresponds to constant A^ larger 
than and thus to scale-independent intermittent distributions as shown in 



lim ^= 

A-*o \ (TAV27r 



exp 



ln'^{a/ao) 
2A2 
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10-1 



10-2 



10-3 



10 



1000 



T [a.u.] 



Figure 4: The form parameter for the data sets OnS, Onl, Off, On2 and Lab are shown in 
double-logarithmic presentation. In that order r = 1 corresponds to 0.25 s, 1 s, 1 s, 0.02 s 
and 10~"'s (0.002 T). The straight line represents a fit according to Eq. 1231 . 

Fig. Inland Fig. [SJa) exemplary for the Onl data set. The A^-values for all data 
sets are illustrated in Fig. IH 

For isotropic turbulence scale-independent distributions occur for large scales 
as well but are always found to be Gaussian [2HI for scales larger than the in- 
tegral time T. From a mathematical point of view these stable Gaussian PDFs 
are the result of a stable stochastic Gaussian process. The other class of stable 
distributions are the so-called Levy distributions - characterised by power-law 
tails - which are the result of a fractional stochastic process. The observed 
atmospheric increment PDFs show robust (stretched) exponential tails that de- 
cay faster than a power-law and slower than a Gaussian distribution. So the 
question arises whether the atmospheric PDFs can be explained as a fractional 
process or as a superposition of different Gaussian processes in analogy to the 
Castaing distribution. To decide this the concept of stability and the connection 
to increment analysis should briefly be introduced. 

3.2.1 Stable Distributions 



Consider the sum = of independent and identical distributed (i.i.d.) 



variables. The variables should be distributed according to a PDF p{x) and the 
PDF of the sum- variable is p{x). The distribution p (or p respectively) is then 
called stable if for large m (m oo) 



is fulfilled jlT] . This means that for sufficiently large m the shape of the distri- 
bution does not change as m increases. 

Transferred to increment analysis an increment over a scale r can be identified 
as variable x and the increment of a larger scale mr as the sum- variable Sm- It 



m 




with X = Ax + B 



(14) 
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is immediately shown that 



m—l 



Umrit) = ^Ur{t + It) 
i=0 



(15) 



which means that a large increment can be expressed as the sum of smaller 
ones. When the PDFs of large and small increments are the same this indicates 
a stable PDF. 

The most famous stable distribution is the Gaussian one. Beside this P. Levy 
j29j showed that there exists a whole class of stable distributions. Restricting 
to symmetric distributions their characteristic functions read: 



For asymmetric distributions the characteristic function becomes more compli- 
cated (e.g. j30jl. The analytical form of the corresponding PDFs is only known 
for some special cases (e.g. for a = 1 it is the Cauchy distribution) but their 
asymptotic behavior is always known and given by: 



This algebraic decay of tails means that all higher order moments larger than 
order a do not exist. Generally, distributions with tails decaying faster than oc 
(defined variance) can only converge to a Gaussian while slower decaying 
tails indicate that they can only converge to a Levy stable law. 

As already mentioned the examined atmospheric PDFs show a faster than 
algebraic decay (compare Fig. |31b) and Fig. |Hla), b), c), d)) so it is expected 
that they converge to Gaussian statistics for large scales. Therefore the ob- 
served robust intermittency should be explained by mixing different Gaussian 
distributions rather than by a fractional stochastic process. 




;0 < a < 2. 



(16) 




(17) 
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3.3 Superposition Model for Atmospheric Turbulence 

In 1^ it was proposed to take the instationarity of long atmospheric time se- 
ries into account. In this sense the observed intermittent form of PDFs for 
all examined scales are found to be the result of mixing statistics belonging to 
different flow situations. These are characterized by different mean velocities 
as schematically illustrated in Fig. 13 When the analysis by means of incre- 
ment statistics is conditioned on periods with constant mean velocities results 
are found to be very similar to those of isotropic turbulence. This can be set 
in analogy to the Castaing distribution that interprets intermittent PDFs as a 
superposition of intervals with different standard deviations. 



u [a.u.] 




■/V>\kHft' 

/Vf^)V^- 



t[a.u.] 

Figure 5: Illustration of different mean velocity intervals. Within these intervals statistics 
should be the same as for isotropic turbulence. The magnitude of variations (standard devia- 
tion) grows with mean velocity according to Eq. 12411 . 

Thus we propose a model that describes the robust intermittent atmospheric 
PDFs as a superposition of those of isotropic turbulent subsets that are denoted 
with p{ur\u) and given by Eq. pij) . Knowing the distribution of the mean 
velocity h{u) the PDFs become: 

oo 

p{Ur) = J du h(u) ■ p{Ur\u) . (18) 


To find a suitable averaging time defining the mean value n is a non-trivial 
problem due to the lack of a distinct mesoscale gap and is not the concern of 
the present paper (see e.g. (ITTj). 

As a first approximation we fix the averaging time to be 10 minutes. In |1] it 
is shown that with this choice a partition into isotropic sequences is achieved. 
Furthermore we assume h{u) to be a Weibull distribution: 




Both assumptions are well established in meteorology [I]. In Fig. El it is shown 
that a Weibull distribution is a good representation of h{u). 
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h(u) 




5 10 15 



u [m/s] 



Figure 6: Symbols represent measured mean velocity distributions (averaged over 10 min) of 
the four atmospheric data sets Onl, On2, On3 and Off. Solid lines are fits according to Eq. 
(IDl . The parameters are: k = 2.0, 7.8, 1.8, 3.3 and A = 2.9, 4.8, 7.1, 8.7 for Onl, On2, 
On3 and Off respectively. 



Inserting Eq. (|19() and Eq. (|11|) into Eq. ()18|) the following expression for 
atmospheric PDFs is obtained: 



k 



Aa2 



oo 




oo 


k / 




























2^2 



-fc-l 



2A2 



(20) 



Parameters j4 and k play a similar role as do and in the Castaing distribution. 



Next we briefly discuss the meaning of the two parameters A and k. The 
former is closely related to the expectation value ((n) = A • r(l + k~^)) of the 
mean velocity. Parameter k determines the form (shape) of the distribution. 
Small values of k correspond to a broad distribution where many different mean 
velocities contribute significantly to the integral in Eq. H20() . Large A:-values 
however prevent this contribution and the distribution is closely centered around 
A. In the limit of very large k a delta distribution is obtained: 



lim 

k—*oo 



k_ 



u 



k-l 



exp 



6iu - A) 



(21) 



So for very large /c-values Eq. (PU]) converges to Eq. ((TT|) i.e. to isotropic 
turbulence. For additionally large r and accordingly vanishing one finally 
gets 



lim lim p{u^ 



1 



cjov 2vr 



exp 



2al 



(22) 
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Figure 7: Standard Deviation as a Function of u and r 



'0 




ojm/s] 



T 



u [m/s] 



a) Standard deviation ctq as a function 
of -It and T - exemplary for the Onl data 
set. 



b) Standard deviation ctq only as 
a function of u (at fixed r) - ex- 
emplary for T — 1 s and r = 50 s. 



which is a pure Gaussian density with variance (Tq(A). For small /c- values the 
distributions remain intermittent even for large scales. Here intermittency is 
caused by the mixing of periods of different mean velocities. 

To apply Eq. 1)20(1 to experimental data one has to know the parameters A, k, 
(To and A'^ . Parameters A and k can directly be estimated from data by fitting a 
Weibull distribution given in Eq. H19|) to the measured distribution of u. This 
is illustrated in Fig. Elfor all examined atmospheric data sets. 
Next = A^(r, It) and = so{t,u) have to be known. Prom isotropic turbu- 
lence it is well-established that the form parameter decreases monotonously in 
scale. Scaling of flatness according to Eq. together with Eq. (|TH|) directly 
leads to a logarithmic depenency of A^ on r: 



To confirm this in Fig. |1] the form parameter of the Lab data set is fitted with 
Eq. (f^ . For a deeper discussion on the scale dependence of A^ see j^H]- 
The dependence of A^ on u is small (e.g. jSj found a slower decay of A^ for in- 
creasing Reynolds numbers Re (x u). As a first approximation we will disregard 
it. 

The parameter (Tq - the most probable standard deviation - generally grows 
with scale and mean velocity. The scale dependence drops out automatically 
when considering normalized PDFs as it is done here. For the u-dependence 
we assume a linear relation 



which is supported from the measured standard deviation fio, shown in Fig. d 
Thus, to apply Eq. (|2U|) we need to calculate A and k from measured mean 
velocity distributions. Additionally we consider a logarithmic dependence of A^ 



= au-bu- ln{T) 



(23) 



So = br ■ U 



(24) 
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on r and a linear one for ctq on u. Given that we are interested in a preferably 
simple reconstruction of PDFs by means of Eq. H2U() we will disregard the more 
complex parameter relations here. 

With these parameters atmospheric increment PDFs can be determined well 
from Eq. H20() . As shown in Fig. El for Onl, On 3 and Off rather small A;- values 
of 2.0, 1.8 and 3.3 are obtained indicating that resulting PDFs are heavy-tailed 
also for large scales. In Fig. andEb) and Fig. |HJi) the corresponding PDFs 
as well as fits according to Eq. (|20j) are presented. The fits agree quite well 
with measured distributions. For all three data sets the tails show a similar 
decay for large scales and are close to a straight line which corresponds to an 
exponential decay due to the semi-logarithmic presentation. 
In contrast to this behaviour the PDFs of On2 and the fits according to Eq. 

- shown in Fig. |Ht) - approach a Gaussian distribution for large scales 
which is in accordance with the large /c-value of 8.7 as it is found in Fig. 1^1 This 
large k-value mirrors the fact that On2 was measured just over a short period 
of time (approximately 1 hr) with rather constant flow conditions. This means 
that as the flow becomes more and more stationary the superposition of differ- 
ent mean velocities becomes weaker and weaker. In the limit of a stationary 
velocity - as realized for Lab and approximately for On2 - isotropic turbulence 
is recovered and large scale intermittency disappears. 

Note, the intention of using Eq. ((^ is not to get the best possible fits of 
the measured PDFs but to model PDFs at specific locations (with a specific 
mean velocity distribution) a priori. The measured PDFs could also be fitted 
by stretched exponential distributions 11 or by Eq. Hll|) as it is done in Fig. 
ISb) for instance. These fits are based on a posteriori measurements however. 
Neither scaling models nor Eq. (jllj) are able to explain a fiatness larger than 
3 and correspondingly intermittent PDFs for large scales outside of the inertial 
range as it is done by Eq. (PU)) . 

This has interesting consequences for atmospheric velocities. For isotropic tur- 
bulence statistical properties of the flow can well be modelled. Thus knowing 
the composition of isotropic subsets it should be possible to model atmospheric 
velocity statistics as well. 
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Figure 8: Atmospheric Probability Density Distributions: Symbols represent the normalized 
PDFs of the four atmospheric data sets in semi-logarithmic presentation. Straight lines corre- 
spond to a fit of distributions according to Eq. I|20^ . All graphs are vertically shifted against 
each other for clarity of presentation. 




a) Ojf: From top to bottom r takes the 
values: 0.2 s, 10 s, 20 s, 200 s and 2000 s. 



b) Onl: From top to bottom r takes 
the values: 0.5 s, 2.5 s, 25 s, 250 s and 
4000 s. 



P(Ut) [a.U.] 




-10 10 



c) On2: From top to bottom r takes 
the values: 2 ms, 20 ms, 200 ms and 
2000 ms. 



p(U'c) [a.U.] 




-10 10 



d) On3: From top to bottom r takes the 
values: 0.25 s, 2.5 s, 25 s and 250 s. 
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4 Conclusions 

The presented analysis provides a good way to estimate and to explain robust 
and markedly intermittent PDFs of atmospheric increments. As desired the 
method recovers results of isotropic turbulence but is also able to take the large 
scale variations into account properly. A rigorous separation of time scales - 
which is very difficult - is not necessary for our model. 

Standard analysis applied to isotropic turbulence cannot be transferred to long 
atmospheric time series unambiguously as it was shown by means of scaling 
behaviour of structure functions and flatness. These methods are restricted to 
the incrtial range where the change of shape can be more or less well reproduced 
and where at the large scale boundary of this range statistics become Gaussian. 
It has been shown that the observed anomalous robust intermittency can be 
explained as a superposition of isotropic subsets in accordance with the concept 
of stable distributions. Intermittency on large scales is found to be the result 
of large scale mixing of isotropic turbulent subsets. 

With our approach intermittent atmospheric PDFs can be approximated for 

any location as long as the mean velocity distribution is known. This is of 
importance for the construction of a wind turbine for instance. For construct- 
ing wind turbines only the turbulence intensity (relation between the average 
standard deviation and the mean velocity) is taken into account which is a 
very inadequate description. The loads are determined by increments and their 
occurrence statistics, therefore a model that reproduces the right increment 
distributions for every location should be preferred. 
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